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* .  INTRODUCTION 

This  quarterly  technical  summary  report  covers  the  work 
performed  during  the  period  April  through  June  1967.  Work  previously 
completed  or  currently  in  progress  is  mentioned  only  as  it  relates 
to  analyses  completed  during  this  reporting  period. 

Analyses  completed,  for  which  results  have  been  reported, 
are  discussed  in  Section  II  under  descriptive  headings.  Section 
III  contains  a  discussion  of  the  support  and  service  tasks  per¬ 
formed  for  in-house  projects  and  for  other  VELA-UNIFORM  participants. 
Appendix  A  is  a  listing  of  those  organizations  receiving  SDL  data 
services  during  this  period. 

II .  WORK  COMPLETED 

A.  Multiple  Coherence  of  Seismic  Noise 

Most  basic  data  processing  techniques  for  signal  enhance¬ 
ment  or  identification  depend  upon  the  structure  of  the  noise  within 
the  seismic  array.  If  some  of  the  coherent  noise  is  due  to  site 
characteristics  such  as  consistently  coherent  noises  from  particular 
directions,  then  techniques  using  multiple  coherence  will  help  tc- 
isolate  these  consistent  linear  relations.  Many  optimum  filters  for 
estimating  the  signal  take  account  of  these  linear  relations  impli- 
citly  *>y  weighting  with  the  inverse  of  the  spectral  noise  matrix. 
However,  one  cannot  tell  whether  the  coherent  noise  involved  is  due 
to  noise  generating  events  which  cannot  be  predicted  or  controlled. 
Thus,  the  filters  must  be  recalculated  over  a  period  of  noise  record¬ 
ing  immediately  prior  to  the  arrival  of  each  single  signal.  Part  of 
the  coherent  noise  generated  within  the  array  may  be  due  to  various 
causal  factors  for  a  particular  array.  If  so,  we  can  learn  something 
about  these  factors  by  examining  the  linear  relations  between  the 
various  array  elements.  A  potential  benefit  here  is  that  a  con¬ 
sistent  linear  model  relating  the  different  sub-elements  would  elim¬ 
inate  the  need  for  computing  a  different  set  of  filter  coefficients  for 
each  event . 


The  multiple  coherence  function  can  indicate  how  many  seis¬ 
mometer  outputs  in  an  array  are  necessary  to  properly  determine  the 
seismic  noise  field.  If  there  are  n  independent  seismic  noise  com¬ 
ponents,  then  the  multiple  coherence  function  would  be  unity  when 
(n  +  1 ) s t  seismometers  are  placed  in  an  array  to  measure  seismic  noise 
records.  If  part  of  the  background  is  composed  of  incoherent  noise, 
then  the  multiple  coherence  function  would  indicate  the  percentage  of 
coherent  noise  present  and  the  number  of  seismometers  necessary  to 
define  this  coherent  noise.  The  filter  relations  determined  by  the 
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multiple  coherence  computations  can  then  be  used  in  array  summation 
to  bring  the  noise  into  destructive  interference. 

This  analysis  does  not  guarantee  that  such  optimum  processing 
is  possible.  For  example,  if  the  noise  and  signal  propagation  charac¬ 
teristics  across  the  array  are  identical,  no  velocity  filtering  scheme 
an  be  expected  to  separate  the  two  even  though  the  multiple  coherence 
might  be  unity. 

The  multiple  coherence  function  is  the  frequency  domain  equi¬ 
valent  of  the  prediction  error  filter  in  time.  If  n  input  seismic 
ices  predict  the  (n  +  l)st  trace  in  an  array  completely,  then  the 
multiple  coherence  will  be  unity  and  a  prediction  error  filter  could 
be  used  to  exactly  predict  this  (n  +  l)st  output.  In  fact,  linear 
filter  relations  derived  by  the  multiple  coherence  program  produce 
an  estimate  of  the  (n  +  l)st  trace  which,  when  subtracted  from  the 
actual  (n  +  l)st  trace,  given  a  prediction  error  trace.  Thus  the 
combination  of  the  filter  derived  in  the  multiple  coherence  program 
and  the  subtraction..operation  produces  a  prediction  error  filter  as 
shown  in  the  following  diagram: 


Prediction  Error  Filter 


This  multiple  coherence  function  was  used  to  estimate  the  degree 
of  predictability  for  the  following: 

-  short-period  noise  field  at  LASA 
long -period  noise  field  at  LASA 

short -period  noise  field  at  vertical  arrays  -  UBSO,  GV-TX 
and  AP-OK 

short-period  noise  field  at  UBSO  and  TFSO 

1.  Short-Period  Noise  Field  at  LASA 

The  multiple  coherences  of  short-period  data  for  both  within 
the  subarray  (intrasubarray)  and  between  subarrays  (intersubarray) 
were  computed,  as  shown  in  Figure  1.  It  was  concluded  from  this 
study  that : 

a.  Intersubarray  noise  at  LASA  shows  a  low  multiple  coherence 
over  all  frequencies,  even  with  9  input  channels. 

b.  Even  the  closest  subarrays  at  LASA  show  little  coherency. 

c.  Intrasubarray  noise  at  LASA  shows  a  high  multiple  coherence 
for  frequencies  of  1  cps  and  lower  when  7  or  more  inputs 
are  used.  However,  the  7th,  8th,  and  9th  inputs  are  only 

h  km  away  from  the  output  trace.  When  the  closest  inputs 
are  2  km  away  from  the  output,  the  multiple  coherences  are 
significantly  less  than  .8  for  all  but  microseismic  fre¬ 
quencies  . 

l 

The  digital  programs  will  show  high  coherence  and  multiple 
\  coherence  when  present,  e.  g.  ,  LONGSHOT  event  has  high  multiple 
\  coherence  over  subarrays.  An  Andreanof  earthquake,  a  weaker 
event  than  LONGSHOT,  shows  multiple  coherences  between  those 
of  LONGSHOT  and  noise  -  only  cases. 

e.  Low  multiple  coherence  is  not  due  to  insufficient  lags. 

The  programs  have  sufficient  resolution  to  handle  noise  pro¬ 
pagation  with  velocity,  1  km/sec  and  less. 

f.  Power  spectra  from  the  seismometers  at  500  feet  agree  with 
those  at  200  feet  only  out  to  1.2  cps.  Peyond  this  frequency 
the  power  spectra  density  of  the  deeper  seismometers  continue 
to  fall,  while  the  power  spectra  densities  from  the  shallow 
seismometers  rise. 

g.  The  multiple  coherence  using  shallow  seismometers  as  inputs 
trying  to  .aatch  the  deeper  seismometer  as  output  falls  rapidly 
beyond  1  2  cps. 

h.  The  intrasubarray  multiple  coherence  using  shallow  seis¬ 
mometers  are  both  inputs  and  output  is  approximately  the  same 
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Figure  1.  Multiple  coherences  versus  frequency  for  both 
inter  and  intrasubarray  noise  recorded  at  LAS A 
on  29  October  1965. 


as  for  intrasubarray  multiple  coherences  which  uses  the 
trace  from  the  500-foot  hole  as  output.  This  result  is 
obtained  in  spite  of  a  striking  difference  in  the  power 
spectra  of  the  shallow  and  deep  seismograms. 

i.  The  intersubarray  multiple  coherence  indicates  that  the 
expected  noise  reduction  from  a  prediction  error  filter 
is  about  1  db  over  the  fitting  interval  and  about  0  db 
outside  the  fitting  interval. 

2.  Long-Period  Noise  Field  at  LASA 

The  multiple  coherence  was  computed  of  long  period  LASA  data 
for  the  vertical  components  only.  These  seismometers  are  located 
at  the  center  of  each  subarray.  Consequently,  all  computations  are 
intersubarray  coherences  with  a  maximum  of  21  channels  available. 
Figure  2  shows  the  multiple  coherences  versus  frequency  for  two  of 
the  three  time  samples  tested.  Results  led  to  the  following  con¬ 
clusions: 


a.  Multiple  coherences  for  all  long-pexiod  samples  increase 
significantly  with  an  increase  in  the  number  of  input 
channels . 

b.  From  mid-range  of  the  long-period  pass  band  toward  the 
microseismic  frequencies  (7-20  seconds  period)  the  multiple 
coherences  are  greater  than  .65  with  8  or  9  input  channels. 

c.  The  expected  noise  reduction  from  a  prediction  error  filter 
in  the  fitting  interval  is  as  much  as  9  db  at  16  seconds 
period.  The  noise  reduction  outside  the  fitting  interval 
at  16  seconds  period  is  from  3-4  db. 

<3.  Over  the  7-20  second  period  range,  the  expected  noise 
reduction  from  a  prediction  error  filter  in  the  fitting 
interval  is  about  4  db.  Over  the  same  7-20  second  period 
range  an  average  of  1  db  is  obtained  outside  the  fitting 
interval . 

e.  The  db  improvement  figures  given  above  were  computed  when 
only  7  input  channels  were  available.  We  estimate  that  the 
db  improvement  both  w_thin  and  outside  the  fitting  interval 
would  be  increased  with  more  input  channels. 

f.  Over  the  7-20  second  period  range  the  multiple  coherences 
for  all  samples  tested  are  quite  similar  to  each  other  up  to 
6  or  7  input  channels.  Since  the  9  channels  available  for 
one  sample  showed  a  significant  increase  in  the  multiple 
coherence  over  the  6  and  7  channel  examples,  we  conclude 
that  at  least  9  input  channels  are  necessary  to  adequately 
model  the  long  period  noise  at  LASA. 
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Figure  2.  Multiple  Coherence  vs.  Frequency  for  two  LP-Z  Samples  from  LASA. 
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Short-Period  Noise  at  Three  Vertical  Arrays  -  UBSO,GV-TX,  AP-OK 

The  multiple  coherence  of  short-period  vertical  component  data 
from  the  three  vertical  arrays  was  computed.  Figure  3  shows  multiple 
coherence  versus  frequency  for  the  two  output  cases.  The  upper 
diagram  shows  multiple  coherence  with  a  surface  instrument  as  output. 
The  lower  diagram  shows  the  multiple  coherence  with  the  deepest  in¬ 
strument  as  output.  It  was  concluded  that: 

a.  The  multiple  coherence  properties  of  the  noise  at  the  three 
vertical  arrays  with  similar  geometries  are  similar.  At 
all  three  sites  the  downhole  channels  correlate  with  each 
other  better  than  with  the  surface  channel.  At  all  three 
sites  the  coherence  between  downhole  seismometers  decreases 
markedly  as  their  separation  increases. 

b  When  the  deep  trace  is  the  output  and  all  seismometers  are 
added  as  inputs,  the  multiple  coherence  is  highest  at  GV-TX. 
In  fact  it  is  greater  than  .85  for  all  frequencies  from  .1 
to  2.5  cps .  For  the  AP-OK  and  UBSO  wells  the  multiple  co¬ 
herence  varies  between  .6  to  .9  for  most  frequencies  greater 
than  .5  cp^. 

c.  All  three  sites  show  a  significant  2  cps  component  in  the 
noise  spectra  for  most  of  the  examples  tested.  The  multiple 
coherence  of  this  2  cps  component  is  high  but  not  appreciably 
higher  than  other  noise  frequencies  above  the  microseismic 
band.  Even  the  surface  trace  correlates  well  with  the  down¬ 
hole  traces  for  the  2  cps  component. 

d.  The  noise  at  AP-OK  appears  to  be  more  stationary  than  noise 
at  the  surface  arrays.  For  three  successive  3*5  minute 
samples  the  expected  noise  reduction  from  a  prediction  error 
filter  applied  in  the  fitting  interval  was  16  to  20  db  over 
the  microseismic  band  and  between  2  and  7  db  for  frequencies 
greater  than  0.5  cps.  Outside  the  fitting  interval  the  ac- 
pected  noise  reduction  was  within  2  db  of  these  values. 

4,  Short-Period  Noise  at  UBSO  and  TFSO 

The  multiple  coherences  were  computed  from  arrays  of  short-period 
vertical  component  seismometers'.  Figures  4  and  5  show  the  multiple 
coherence-  versus  frequency  for  two  noise  samples  recorded  at  UBSO 
and  TFSO  respectively.  It  was  concluded  that: 

a.  The  multiple  coherence  of  the  noise  at  UBSO  and  TFSO  short- 
period,  vertical  component  arrays  is  high  (greater  than  0.9) 
over  the  microseismic  frequency. 

b.  The  decay  of  the  multiple  coherence  of  the  noise  with  in¬ 
creasing  frequency  is  faster  at  TFSO  than  at  UBSO  and 
faster  at  UBSO  than  at  LASA  subarrays . 
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Figure  3.  Multiple  coherence  versus  frequency  for  a  UBSO 

sample  recorded  in  December  1966.  The  upper  diagram 
uses  the  surface  instrument  as  output;  the  lower 
uses  the  deepest  instrument  as  output. 
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Figure  5.  Multiple  Coherences  vs.  Frequency  with  2-9  Inputs 

for  two  noise  samples  recorded  at  TFSO,  January  1967 
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The  noise  has  a  low  multiple  coherence  (less  than  .4) 
with  up  to  9  inputs  at  UBSO  and  TFSO  for  frequencies 
greater  than  1.5  cps .  Thus  the  noise  is  largely  incoherent 
over  the  upper  part  of  the  signal  band  at  UBSO  and  TFSO. 

d.  The  noise  at  UBSO  and  TFSO  over  adjacent  time  samples 

is  stationary  only  over  the  microseismic  frequency  band. 

B,  LASA  Analysis 

1.  The  Detection  Threshold  at  The  Montana  LASA 

As  part  of  the  continuing  stud^  to  evaluate  the  Montana  LASA,  an 
analysis  was  made  to  determine  a  detection  threshold  at  the  LASA  for 
events  originating  at  teleseismic  distances.  Any  such  magnitude  value 
is  highly  generalized  since  it  depends  not  only  cn  the  epicentral  dis¬ 
tance  and  azimuth,  depth  of  the  hypo-center,  and  the  source  mechanism, 
but  also  on  the  precision  of  the  beamforming  process  and  the  type  of 
prefiltering  used.  Moreover,  it  is  undoubtedly  a  function  of  the  re¬ 
cording  period.  By  this  is  meant  that  the  nighttime  noise  level  is 
lower  than  that  recorded  during  the  daylight  hours,  and  detection 
threshold  estimates  will  be  lower  than  that  obtained  from  a  similar 
set  of  events  recorded  during  the  day.  For  this  study,  the  basic  pro¬ 
cedure  was  to  use  LASA  beamfo .med  outputs  from  56  events  recorded  at 
night  to  compute  event  magnitudes  and  corresponding  threshold  values. 

Table  1  summarizes  the  results  of  this  analysis  and  shows  C&GS 
magnitude,  LASA  Magnitude,  LASA  threshold  magnitude,  and  averages  for 
the  56  events  analyzed  to  date.  These  data  illustrate  two  important 
points.  First,  the  average  C&GS  magnitude  differs  from  the  average 
LASA  magnitude  by  less  than  0.1,  even  though  LASA  computations  were 
based  on  data  taken  from  beamformed  output  traces.  Second,  the 
Montana  LASA  should  be  able  to  detect  a  mixture  of  shallow  and  deep 
events  of  average  magnitude  3.6  as  a  result  of  the  beamforming  process. 
For  bombs  and  shallow  quakes,  the  threshold  would  be  0.1  higher  be¬ 
cause  of  the  correction  factor  for  deep  foci.  For  daytime  recordings 
we  estimate  that  it  would  be  an  additional  0.1-0. 2  higher  because  of 
higher  noise  levels.  Obviously,  this  magnitude,  3.6,  is  highly  general¬ 
ized  since  it  is  a  function  of  several  parameters,  which  include  epi¬ 
central  distance  and  azimuth,  source  mechanism,  beamforming  precision, 
and  type  of  bandpass  filter  used. 

2.  Distribution  of  Short-Period  P- Phase  Amplitudes  Over  LASA 

In  the  following,  the  distribution  of  the  peak  to  peak  (p-p) 
amplitudes  for  several  events  at  LASA  will  be  considered  and  it  will 
be  shown  that  this  data  approaches  a  log-normal  distribution. 

Plotting  the  p-p  amplitudes  of  all  SPZ  seismometers  at  LASA  for 
a  single  event  against  frequency  of  occurance  yields  a  graph  which  seems 
to  have  no  familiar  probability  distribution.  Figure  6,  which  shows 
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Table  1.  Event  Magnitudes  at  LASA 


the  p-p  amplitudes  of  an  earthquake  is  a  typical  graph  of  this  type. 
However,  if  the  p-p  amplitudes  that  are  associated  with  a  given 
subarray  are  inspected,  a  pattern  emerges.  There  is  generally  a 
modal  clustering  with  a  short  tail  toward  decreasing  values  and  a 
long  tail  toward  increasing  values  suggesting  that  the  data  may  come 
from  a  log-normal  distribution  (Figure  7). 

It  was  evident  that  some  subarrays  exhibit  a  higher  degree  of 
sensitivity  for  events  from  one  source  region  than  for  events  from 
another.  These  consistencies  indicate  that  the  individual  subarrays 
are  biased  in  a  certain  sense  due  to  the  travel  path  of  the  source 
signal.  Since  the  elements  in  each  subarray  are  closely  spaced, 
in  comparison  to  the  distance  separating  subarray  centers,  it  was 
found  that  the  individual  subarray  preference  from  all  of  the  subarrays 
could  be  eliminated  by  considering  each  subarray  as  a  sampling  of  the 
population  of  the  p-p  amplitudes  and  standardizing  each  sampling  in 
the  usual  sense,  i.e.,  if  a  subarray  has  a  N(m,s)  distribution  des¬ 
cribing  the  behavior  of  the  observations  x,  then  the  linear  transfor¬ 
mation  will  give  the  standardized  distribution  N(0,1). 


By  applying  this  transformation  to  all  subarrays  over  every  event,  all 
subarrays  and  all  events  can  be  pooled  and  then  the  properties  of 
the  entire  pool  can  be  investigated.  If  the  assumption  that  the  log¬ 
arithms  of  the  p-p  amplitudes  in  each  subarray  are  normally  distri¬ 
buted  is  valid,  then  the  pooled  distribution  will  also  be  normally 
distributed . 

One  hundred  and  fifty  seven  seismograms  of  eight  events  yielded 
3,925  usable  amplitude  records  which  were  used  to  test  this  hypothesis. 
Results  showed  a  definite  tendency  toward  normality  of  the  data.  Plots 
of  the  data  gave  a  visual  indication  regarding  the  normality  of  the 
data.  Statistically,  the  distribution  was  \2  -tested  for  normality 
with  the  result  that  at  the  5%  level  of  significance  the  sample  dis¬ 
tribution  is  consistent,  with  the  hypothesis  that  the  parent  distri¬ 
bution  is  normal.  • 

3.  Correlation  of  Amplitude  Anomalies  at  LASA 

Certain  consistencies  have  been  shown  to  exist  among  the  anomalies 
of  short  period  teleseismic  P-phase  amplitudes  when  those  anomalies 
are  grouped  according  to  geographical  regions.  The  purpose  of  this 
study  was  to  present  evidence  that  the  source  path  strongly  influences 
the  observed  amplitude  anomalies. 

Eight  earthquakes  from  the  same  geographical  area,  viz.,  the  Fiji 
Islands,  were  selected  for  this  investigation.  All  eight  events 


Figure  6.  P-P  Amplitudes  of  an  Earthquake 
Recorded  at  All  LASA  SPZ  Seismoi 


occurred  at  143°  Az  with  respect  to  the  center  of  the  LASA  and  ranged 
from  9,500  km.  to  10,500  km.  distance.  These  events  were  ordered 
in  relation  to  their  distance  from  the  LASA  and  were  partitioned  into 
pairs.  The  result  of  this  grouping  was  a  set  of  four  pairs  of  events 
such  that  the  elements  within  each  pair  had  epicenters  within  eight 
kilometers  of  each  other.  Although  these  epicenters  were  close  to¬ 
gether,  the  focuses  of  the  events  had  a  considerably  greater  range. 

Each  of  the  pairs  was  investigated  for  the  presence  of  correlation 
between  member  events  in  the  following  manner.  The  peak  to  peak 
P-phases  amplitudes  recorded  by  the  center  elements  of  the  subarrays 
were  determined  for  each  event.  In  this  way,  an  ordered  set  of  ampli¬ 
tudes  indexed  by  the  recording  instruments  was  associated  with  each 
event.  Thus,  in  order  to  detect  the  presence  of  correlation  between 
any  two  events,  one  merely  tests  for  correlation  between  the  two 
sets  of  amplitudes  generated  by  respective  events.  But  before  taking 
this  step,  it  is  important  to  note  that  an  accurate  interpretation 
of  the  correlation  coefficient  is,  in  general,  possible  only  when 
the  underlying  population  is  normally  distributed.  It  was  demonstrated 
in  the  discussion  under  2,  above,  that  a  log-normal  distribution  des¬ 
cribes  the  behavior  of  amplitude  anomalies  over  LASA  with  remarkable 
precision.  Hence,  by  using  the  logarithms  of  the  peak  to  peak  ampli¬ 
tudes,  each  event  has  been  characterized  with  a  normally  distributed 
set  of  measurements  and  the  correlation  theory  can  now  be  applied. 

The  estimate  of  the  coefficient  of  correlation  is  defined  by 

Cov  (X ,  Y) 

r  =  - -  -  . 

7~Var  X  "  Var  Y 

Results  showed  that  one  pair  was  correlated,  one  was  borderline, 
and  two  were  uncorrelated.  The  pairs  were  ranked  from  the  highest 
to  the  lowest  degree  of  correlation;  then  the  pairs  were  ranked  from 
those  with  earthquakes  which  occurred  closest  together  to  those  with 
earthquakes  farthest  apart.  These  two  rankings  yielded  the  same  order 
in  both  cases.  The  events  which  occurred  closest  together  had  the 
highest  correlation  and  those  which  occurred  farthest  apart  were  least 
correlated . 

Similar  tests  were  performed  on  three  teleseismic  events  which 
occurred  near  the  surface,  and  at  the  same  location.  All  possible 
combinations  of  pairs  were  tested  and  in  each  instance  a  high  positive 
correlation  was  found  to  exist. 


The  underground  nuclear  explosion  BILBY  was  detonated  in  the 
Nevada  Test  site  on  13  September  1963.  BILBY  generated  seismic  waves 
which  were  recorded  well  at  most  of  the  North  American  stations  and 
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especially  at  the  LRSM  stations.  Furthermore,  it  was  followed  by  the 
cavity  collapse  (assigned  magnitude  =  4.5).  The  presence  of  this 
collapse,  the  relatively  large  size  of  the  explosion,  (m  =  5„8) , 
radiation  of  both  Rayleigh  and  Love  waves,  and  the  proximity  of  the 
site  to  those  of  other  well  recorded  explosions  such  as  HAYMAKER 
and  SEDAN  have  made  BILBY  ideally  suited  for  a  study  to  determine  the 
mechanism  of  generation  of  seismic  waves  and  the  radiation  pattern. 

Because  of  the  particular  nature  of  the  desired  study  and  to 
insure  that  a  procedure  of  analysis  similar  to  that  previously  used 
was  exactly  followed,  it  was  decided  to  have  Dr.  M.  Nafi  Toksoz, 
Massachusetts  Institute  of  Technology,  conduct  this  analysis,  assist- 
ed  by  Mr.  Kevin  Clermont,  Princeton  University. 


Basically,  the  procedures  followed  were  to  first  minimize  the 
effects  of  the  propagation  by  determining  the  radiation  patterns 
of  the  P-waves  and  Rayleigh  waves  of  BILBY  relative  to  those  of  the 
collapse.  Then  the  ratios  of  the  Love  and  Rayleigh  wave  amplitudes 
were  used  to  determine  a  source  mechanism  which  was  consistent  with 
the  Rayleigh  wave  radiation  pattern.  The  source  time  function  was 
determined  from  the  amplitude  spectra  corrected  for  propagation  and 
instrument  effects.  Finally,  the  BILBY  results  were  compared  with 
the  HARDHAT ,  HAYMAKER,  SEDAN,  and  SHOAL  explosions.  The  most  signi¬ 
ficant  result  is  that  BILBY,  like  the  HAYMAKER  and  SHOAL  explosions, 
generated  some  Love  waves.  The  source  mechanism  in  all  cases  can 
be  explained  in  terms  of  an  ideal  explosive  source  superimposed  over 
a  tectonic  source  of  double-couple  form.  The  orientation  and  relative 
strength  of  the  double-couple  seem  to  be  controlled  by  the  properties 
of  the  medium  and  the  orientation  of  the  tectonic  axis.  BILBY  (in 
tuff)  aii<^  HAYMAKER  (in  alluvium)  were  located  about  5  km  apart.  The 
radiation  patterns  of  Rayleigh  waves  are  almost  identical.  In  both 
cases  the  principal  plane  of  the  double-couple  is  oriented  in  the 
dilection  Q  -  340°.  For  the  SHOAL  explosion,  which  was  detonated 
in  a  completely  different  area,  the  orientation  of  the  double-couple 
was  ^-n  verY  agreement  with  those  of  earthquakes  in  the  area. 

These  facts  suggest  that  multipolar  contributions  to  the  radiation 
patterns  are  controlled  by  the  general  tectonic  features  of  the 
region . 

D .  Synthesis  of  Additive  Ambient  Seismic  Noise  With  A  Gaussian 
Markov  Model  '  “ 

The  purpose  of  this  study  was  to  show  that  simple  synthesized 
seismic  noise  models  can  qualitatively  account  for  many  of  the  observed 
properties  of  the  noise.  Such  models  are  useful  for  efficiently 
generating  standard  test  data  on  a  digital  computer  for  checkout  of 


vertical  array  processors  and  for  preliminary  evaluation  and  ranking 
of  seismic  signal  processors. 

The  ambient  seismic  noise  was  modeled  by  a  single  Gaussian 
population  from  whicii  independent  realizations  or  states  were  taken 
as  input  to  tuned  filters  with  spectral  peaks  matched  to  those  ob¬ 
served  in  noise  samples,  for  example  at  .2  cps  and  2.  cps.  For 
each  spectral  noise  peak,  the  realization  on  channel  i  +  1  was 
equal  to  a  constant  times  the  realization  on  channel  i  plus  another 
constant  times  a  new  realization  on  channel  i  +  1.  Results  showed 
that  the  constants  defining  the  Markov  process  could  be  used  to 
theoretically  derive  the  associated  power  spectral  matrix  of  the 
noise  model.  The  model  could  be  extended  to  dispersive  systems  by 
using  a  set  of  constants  and  time  lags  to  relate  the  noise  on 
channel  i  to  that  on  channel  i  +  1. 

A  vertical  array  signal  model  was  also  given.  The  purpose  was  to 
efficiently  generate  noise  and/or  signals  at  prescribed  S/N  ratios. 

The  noise  covariance  structure  was  close  to  that  observed  naturally 
and  was  known  exactly  for  the  noise  model  realizations.  Thus  the 
spectral  covariance  of  the  noise  was  given  exactly,  subject  only  to 
roundoff  error,  and  conditions  of  stationarity  and  equilibrium  were 
satisfied  by  the  data  generated  for  testing  and  designing  multi¬ 
channel  filters. 

E .  Principal  Component  Analysis  of  Seismic  Data 

Two  technical  notes  were  sxibmitted  to  the  SDL  by  Measurement 
Analysis  Corporation  as  part  of  specialized  studies  which  they  per¬ 
formed  for  the  SDL.  The  work  was  concerned  with  the  development  of 
the  principal  component  theory  to  seismic  noise  records. 

Principal  component  analysis  is  designed  to  explain  observed 
relations  among  sample  records  obtained  from  an  array  of  seismometers 
in  tt  ms  of  simpler  relations  by  reducing  the  random  variables  of 
these  records  to  a  small  number  of  linear  combinations.  The  sum  of 
the  variances  of  all  principal  components  is  the  sum  of  the  variances 
of  the  original  variables.  Thus,  if  there  exist  principal  components 
with  large  variances  which  account  for  most  of  the  variability,  the 
dimensionality  of  the  problem  might  be  reduced  by  attention  only  to 
these  principal  components. 

The  first  note  concerned  itself  with  the  development  of  the 
principal  component  theory  by  first  summarizing  a  statistical  inter¬ 
pretation  of  the  principal  component  and  then  by  analysis,  based  on 
noise  data  recorded  at  the  Montana  LASA  The  analysis  consisted  of 
obtaining  the  proportion  of  the  total  variance  (power)  from  10  seis¬ 
mometers  as  explained  by  the  first  four  principal  components.  Results 


showed  that  at  0.2  cps  approximately  80%  of  the  total  variance  was 
accounted  for  by  the  first  principal  component. 

The  second  note  concerned  itself  with  utilizing  the  principal 
component  to  determine  the  noise  source,  based  on  data  from  LONGSHOT 
and  a  geographically  nearby  ef^thquake  recorded  at  a  LASA  subarray. 
The  phase  shifts  of  the  first  principal  component  at  0.2  cps  cor¬ 
rected  for  LASA  instrument  response  resulted  in  an  accurate  estimate 
of  the  general  direction  of  the  main  noise  source.  Computed  examples 
of  the  data  verified  this  estimate. 


Measurement  Analysis  Corporation  presented  a  report  which  des¬ 
cribed  some  interpretations  and  uses  of  eigenvalues  and  eigenvectors 
of  spectral  and  sample  spectral  density  matrices  of  multiple  stationary 
time  series. 

The  spectral  density  matrix  of  a  zero-mean  multiple  stationary 
time  series  was  defined.  Eigenvalues  and  eigenvectors  of  the  spectral 
density  matrix  were  discussed  and  principal  component  theory  was 
presented.  Statistical  distribution  theory  and  related  results  were 
used  to  investigate  the  eigenvalues  of  a  sample  spectral  density 
matrix.  This  investigation  gave  methods  for  obtaining  simultaneous 
confidence  bounds  on  the  elements  of  the  true  spectral  density  matrix 
and  its  inverse,  and  also  methods  for  obtaining  confidence  bounds 
on  the  eigenvalues  of  the  true  spectral  density  matrix. 

G .  Shot  and  Earthquake  Analyses 

Standard  statistical  summary  reports  were  issued  on  the  Nevada 
Test  Site  nuclear  t vents  BOURBON,  NASH,  and  GREELEY.  Table  2  sum¬ 
marizes  part  of  the  seismic  data  compiled  from  these  analyses. 

Other  projects  completed  deluded  photo  albums  from  four  events 
and  travel-time  computations  for  NTS  and  other  sites. 


HI.  SUPPORT  AND  SERVICE  TASKS 

A .  VELA-UNIFORM  Data  Services 

As  part  of  the  contract  work-statement,  the  SDL  pro¬ 
vided  one  or  more  of  the  following  support  and  service  functions 
for  VSC  and  other  VELA  participants: 

-  copies  of  16  and  35  mm  film 

-  playouts  of  earthquakes  and  special  events 
copies  of  existing  composite  analog  tapes 

-  composite  analog  tapes  of  special  events 

-  use  of  1604  computer  for  checking  out  new  programs 
or  running  production  programs 

-  copies  of  digital  programs 

-  digitized  data  in  standard  formats  or  special  formats 
for  use  on  computers  other  than  the  1604 

-  running  SDL  production  programs,  such  as  power  spectral 
density  and  array  processing  cm  specified  data 

-  digital  x-y  plots  of  power  spectra  or  digitized  data 
signal  reproduction  booklets 

-  copies  of  MIT  Geophysics  Program  Set  II 

space  for  visiting  scientists  utilizing  SDL  facilities 
to  study  data  and  exchange  information  with  SDL  person¬ 
nel  .  . 

t.  . 

During  this  report  period,  51  such  projects  were  com; 
pleted  and  the  15  organizations  receiving  these  services  are  listed 
in  Appendix  A. 


B.  Data  Library 

The  Data  Library  contains  approximately  7,000  digitized 
seismograms,  170  digital  computer  programs  and  290  composite  analog 
magnetic  tapes,  all  available  for  use  by  the  VELA-UNIFORM  program. 

The  following  additions  were  made  during  this  period. 

1.  Digital  Seismograms  -  372  including 

-  data  from  10  explosions  and  3  underwater  events 

-  noise  samples  from'LASA,  TFO,  UBSO,  CPSO,  and  WMSO 
noise  from  several  vertical  arrays 

-  earthquake  data  recorded  at  various  stations 

2.  LASA  Data  -  84  digital  tapes 

-  there  are  a  total  of  785  digital  tapes  in  the  library. 

There  is  also  a  master  calibration  tape  which  contains 
the  magnification  (digital  counts  per  millimicron)  of 
each  sensor  ‘for  every  subarray.  These  magnifications  have 


been  computed  for  all  calibration  tapes  currently  in  house. 

As  each  new  calibration  is  received,  it  is  routinely  run 

through  the  program  CALI5R  and  added  to  the  master  tape . 

3.  Digital  Programs  -  19  including: 

ARCTANGENT  FUNCTION  -  to  computa  the  arct agent  of  the  ratio 
X/Y,  based  on  certain  conditions. 

CALIBR  —  to  compute  the  magnification  level  (counts  per 
millimicron)  of  each  LASA  seismometer 

PARDATA  -  to  take  an  input  tape  containing  a  label  and  i 
no.  of  data  traces  X^ ,  partition  X^  into  j  no.  of  data  traces 

X^,  and  write  a  new  tape  containing  either  a/or  no.  label  and 
il  i2  13 

RELABLIB  -  using  standard  FORTRAN-63  library  seismogram 
tapes,  either  the  packed  or  unpacked  format,  this  program  will  relabel 
seismograms,  copy  seismograms,  skip  seismograms,  merge  tapes,  or  any 
combination  of  these  operations. 

\  v  .  \ 

SPEC MAT  -  to  compute  the  theoretical  power  spectral  density 
matrix  for  the  Gaussian  Markov  model. 

SPECMTRX  -  to  compute  theoretical  power  spectral  densities, 
coherencies,  and  phases  for  the  Gaussian  Markov  chain. 

LOPER  -  to  reformat  the  long  period  multiplexed  digital  LASA 
data  which  were  recorded  at  5  points  per  second  on  the  slow  mode 
tapes  into  standard  SDL  library  tape  format. 

SHOPER  -  to  reformat  the  short  period  multiplexed  digital  LASA 
data  which  were  recorded  at  20  points  pe_  second  on  the  slow  mode  ' 
tapes  into  standard  SDL  library  format . 

FORDEC  -  to  read  a  standard  SDL  library  tape  whose  sampling 
rate  is  20  samples  per  second  and  decimate  the  data  to  10  samples  per 
second  and  write  it  onto  a  tape  in  standard  SDL  library  format  with 
corrections  in  the  label  of  the  new  tape  for  sampling  rate  and 
number  of  points. 

RMS  ONLY  -  to  read  up  to  9  seismograms  from  a  magnetic  tape, 
in  subset  tape  format,  detrend,  filter  (optional),  compute  rms  on  a 
particular  desired  region,  and  the  mean  of  the  rms. 


LOPDEC  -  to  reformat  the  long  period  multiplexed  digital  LASA 
tapes  which  were  recorded  at  5  points  per  second  into  standard  SDL 
library  tapes  decimated  to  1  point  per  second. 


SPEC SUB  -  to  read  up  to  6000  points  of  data  in  sub-set  format, 
detrend,  demagnify,  filter  it  and  form  a  new  sub-set  tape. 


LAS AN  -  to  read  fast  mode  LASA  data,  20  samples  per  second, 
less  than  3000  points,  from  a  standard  SDL  library  tape  and  further 
process  the  data. 

ALAST  -  ALAST  reads  the  save  tapes  prepared  by  MODLAS  or 
LASAN  and  aligns  the  traces  according  to  the  Delta  T's  and  Anomalies 
supplied  by  the  user  and  makes  a  new  save  tape  with  traces  aligned 
and  in  sub- set  format . 

MODLAS  -  to  read  LASA  data  from  an  SDL  library  tape  and  further 
process  data. 

LOPSAN  -  LOPSAN  is  a  1604  program  designed  to  read  LP  data 
from  SDL  library  format  tapes  whose  sampling  rates  are  1  to  5  samples 
per  second,  up  to  12000  points,  and  further  process  the  data. 

PREDICT  -  to  compute  a  prediction  of  the  straight  sum  of  Multi¬ 
channel  arrays,  ispan  units  ahead  in  time  where  ispan  is  the  predict¬ 
ion  span .  The  predicted  sum  and  the  actual  sum  are  subtracted  to 
yield  a  prediction  error,  and  all  are  plotted  to  the  same  scale  factor. 
Also,  to  compute  and  plot  power-spectra  of  these  traces. 

VFKSPTRM  -  to  compute  and  display  the  frequency-wave  number 
power  spectra  of  seismic  noise  along  with  a  response  function  for  the 
corresponding  vertical  array. 

AMBILOG  -  to  convert  a  Geotech  Ambilog  digital  tape  into  the 
standard  SDL  unpacked  library  tape  format  and  produce  a  seismogram 
label  listing  and  punched  index  cards. 

4.  Analog  Composite  Tapes  -  10  including: 

a.  Made  by  SDL 

-  OONW  data  for  Lincoln  Laboratory 

-  Earthquake  data  for  USCGS 

b.  Made  by  Geotech 

-  PERSIMMON 

-  AGILE 

-  14  - 
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C .  Data  Compression 

This  is  a  continuing  routine  operation,  and  production  is 
maintained  at  the  level  needed  to  meet  the  requirements  of  the 
field  operation  (LRSM  and  U.  S.  Observatories)  and  the  Seismic  Data 
Laboratory.  For  this  period,  1,244  tapes  were  compressed. 


D .  Automated  Bulletin  Process 

January,  February,  and  March  1967  LRSM  and  Observatory 
bulletins  were  processed  during  this  report  period  and  forwarded  tc 
the  Geotech,  A  Teledyne  Company,  for  checking  and  publication. 


APPENDIX  A 


ORGANIZATIONS  RECEIVING  SDL  DATA  SERVICES 
April  -  June  1967 

Colorado  School  of  Mines 
Earth  Sciences,  Teledyne 
General  Atronics  Corporation 
Geotech,  Teledyne 
IBM  Corporation 

Lamont  Geophysical  Observatory 
Massachusetts  Institute  of  Technology 
Oregon  State  University 
Pennsylvania  State  University 
Princeton  University 
San  Calixto,  Bolivia,  Observatory 
Texas  Instruments,  Incorporated 
University  of  California 
University  of  Michigan 


U.  S.  Coast  and  Geodetic  Survey 


_  Security  Classification 


„  DOCUMENT  CONTROL  DATA  •  R&D 

- l9»Oirt,y  cloooituoUon  o I  My  „/  .b.rr.c,  *,d  mu„  6.  ,h,  ,1DO„ 

I  -  ORIGINATING  ACTIVITY  (Corpotmt*  mulhor)  ~  _ is.  .  . 

^  *•  report  SECURITY  C  LAtSIFtCA  TION 

TELEDYNE ,  INC .  Unclassified _ 

ALEXANDRIA ,  VIRGINIA  22  314  **  G"°U* _ 

*•  report  title  ——  - - - 


[  SEISMIC  DATA  LABORATORY  QUARTERLY  TECHNICAL  SUMMARY  REPORT 

«•  DESCRIPTIVE  NOTES  (Typo  ot  toptrt  and  Inchttlrt  dot*,)  .  . .  . 

Quarterly  Summary  -  April  -  June  1967 

S'  AUTHOUfJ)  (Loot  nomt,  Hnt««u,  I nltlal)  "  — — ...  . .  .  . 

Dean,  William  C. 


report  date 


To-  TOTAL  NO.  OP  PAflKI 
26 


M-  ORIOINATOR'i  REPORT  NUMRCRfJJ 


17  Julv  1967 


•  CONTRACT  OR  SRANT  NO. 

F  33657-67-C-1313 

S.  PROJECT  NO. 

VELA  T/6702 

ARPA  Order  624  g/T,M“J2(,y*0*T 

4.  ARPA  Program  Code  No.  5810  Technical  Summary  Report  No.  16 

To.  AVAIL  ABILITY/LIMIT  A  TIOft  NOTICES  “™ - - - - — - - — 

This  document  is  subject  to  special  export  controls  and  each  trans¬ 
mittal  to  foreign  governments  or  foreign  national  may  be  made  only 
with  prior  approval  of  Chief.  AFTAC . 


II.  SUPPLEMENTARY  NOTES  12.  SPONSORING  MILITARY  ACTIVITY  ——————— 

-  ADVANCED  RESEARCH  PROJECTS  AGENCY 

NUCLEAR  TEST  DETECTION  OFFICE 
- - -  WASHINGTON .  D .  C  .  ' 

II.  ABSTRACT  - - * - — — - - - - 

This  report  discusses  the  work  performed  by  SDL  for  the 
period  April  through  June  1967,  and  is  primarily  concerned 
with  seismic  research  activities  leading  to  the  detection  and 
identification  of  nuclear  explosions  as  distinguished  from 
earthquake  phenomenon.  Also  discussed  are  the  data  services 
performed  for  other  participants  in  the  VELA-UNIFORM  project. 


DD  .Wi.  1473 


issified _ 

Security  Classification 


mmmM 


- Unci ass Lfl ed 

Security  Classification 


14 

KEY  WORDS 

L  INn  A 

.NK  li 

I.IUK  C  1 

’■OLL 

WT 

»OI.E  |  *1 

P-OUE 

*  r 

Seismic  Data  Laboratory  -  Quarterly 

Technical  Summary 

VEI A-UNIFORM  Project 

Seismic  Data  Analysis 

Large  Aperture  Seismic  Array 

Seismic  Noise 

Principal  Components  of  Seismic  Noise 

Coherency  Analyses 

INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee.  Department  of  De¬ 
fense  activity  or  other  organization  (corporate  author)  issuing 
the  report. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether  . 
“Restricted  Data"  is  included.  Marking  is  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

26.  GROUP:  Automatic  downgrading  is  specified  in  DoD  Di¬ 
rective  5200. 10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ised. 

3.  RETORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  caonot  be  selected  without  classifica¬ 
tion,  show  title  classi.icatiOn  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  name(s)  of  authoKs)  as  shown  on 
or  in  the  report.  Entei  last  name,  first  name,  middle  initial. 

If  military,  show  rank  end  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  yean  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  be.,  enter  the 
number  of  pages  containing  information. 

7b.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  *>propriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  waa  written. 

8b,  8c,  8s  8<f.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR’S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

9b.  OTHER  REPORT  NUMBERS):  If  the  report  has  been 
assigned  any  otheT  report  numbers  (either  by  the  originator 
or  by  the  aponaor),  ilso  enter  this  number! s). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  or  this 
report  from  DDC” 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  OtheT  qualified  DDC 
users  shall  request  through 

If 

-  a 

(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC  Other  qualified  users 
shall  request  through 

II 

(5)  "All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 

•  I 

— —  ■  _ -  * 

If  the  report  has  been  furnished  tc  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

1L  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay¬ 
ing  for)  the  reeesrch  and  development.  Include  address. 

13-  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  id  the  body  of  the  technical  re¬ 
port.  If  additional  apace  is  required,  a  continuation  sheet  shall 
be  attached. 

it  is  highly  de'I-sble  thst  the  abstract  of  classified  reports 
be  unclassified  Each  paragraph  of  the  abstract  shall  and  with 
an  indication  i  the  military  security  classification  of  the  in¬ 
formation  in  tl  paragraph,  represented  as  (TS),  (S).  (C),  or  (V). 

There  is  no  limitation  cn  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional 


gpo  see- ssi 

Unclassified _ 

Security  Classification 


